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We report about a simulation using three dynamical Wilson quarks and on the progress in going to small quark 
masses. 



1. INTRODUCTION 

One of the main purposes of lattice QCD sim- 
ulations is to predict the spectrum of QCD, i.e. 
masses and decay constants of the light hadrons. 
Due to computational and algorithmic limitations 
certain simplifications have to be made. The 
quenched approximation allows to reproduce the 
physical spectrum already up to systematic errors 
of 10%, as it was shown in detail in 0. 

Lots of effort has been spent to increase the 
available computing power, see e.g. for such a 
project, and to overcome algorithmic limitations 
js). There have been several large scale simu- 
lations of QCD with two dynamical, degenerate 
quarks Q, all based on the Hybrid Monte Carlo 
algorithm ||^ . Still these simulations worked with 
relatively large quark masses. For phenomenolog- 
ical reasons it is crucial to go beyond this scenario 
and simulate with three dynamical quarks with 
masses eventually approaching the physical ones 

§. _ 

Some restrictions can be overcome by using 

variants of the multibosonic algorithm like 
the two-step multibosonic (TSMB) algorithm 
The latter has already been used successfully in 
simulations of supersymmetric field theories Q 
and finite density QCD g. 

In this contribution we present preliminary re- 
sults from a first simulation using three dynami- 
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cal Wilson quarks. The aim of these simulations 
is to achieve small quark masses and to estimate 
their infiuence on the spectrum. 



2. SMALL MASSES 

There are several problems when tuning the al- 
gorithms to small quark masses. One problem 
is the critical slowing down, which reflects the 
decreasing efficency at smaller masses. In jioj 
it has been shown that this problem is less se- 
vere with multibosonic algorithms than with the 
HMC algorithm. But still this can become a prob- 
lem at too small masses. A first orientation can 
be obtained by making a comparison with pre- 
vious simulations. In jl^ TSMB has been used 
for the = 1 supersymmetric Yang-Mills the- 
ory where the number of flavour is Nf — \- In 
that simulation condition numbers of O(IO^) oc- 
curred. Because the performance of the algo- 
rithm depends mainly on the condition number 
A/e of the squared fermion matrix, and not on 
Nf^ one can similarly expect to reach condition 
numbers of at least 0(10^) for the case Nf = 3. 
This would correspond to quark masses of roughly 
\ms. The polynomial orders ni and n2 that are 
needed for the TSMB algorithm mainly depend 
on A/e. Changing from Nf = ^ to Nf = 3, rii 
is increasing by less than 50% and 712 by about 
25%. This dependence is reflected by the asymp- 
totic estimate of n2 for large A/e |12 
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where Cn; changes only sUghtly with Nf and toq 
is the bare quark mass. 

Another potential problem when approaching 
the chiral limit may be posed by a sign change 
of the fermionic determinant. This so-called sign 
problem may spoil the statistical signal of ob- 
servables. Experience shows that both situa- 
tions (with and without sign problem) are pos- 
sible jljjl^. The sign of the determinant can 
be controlled by the spectral flow of small eigen- 
values. From this experience it is already clear 
that a sign change can only occur if there are ex- 
tremly small eigenvalues (0(10^*)). Such eigen- 
values did not appear in our runs so far, which 
gives a first evidence for the absence of the sign 
problem. 

3. PARAMETER SPACE 

The spectrum of the light hadrons can be stud- 
ied in the confined phase, the phase below the 
finite temperature phase transition 
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lattice: 8x8x8x4 
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Figure 1. Finite temperature phase transition 



On small lattices the transition between the two 
phases can be seen clearly, as e.g. in figure g. The 
parameters of our simulations lead to 



where ro is the Sommer scale parameter. There- 
fore we are safe from finite size effects due to de- 
confinement by choosing lattice sizes L such that 

L/ro > 4. (4) 

The parameters we finally want to achieve 
are shown in table |l|. These parameters give a 
constant volume L/rp — 4 and are such that 
{r^mpsY — 0.8, corresponding to a quark mass 
of TUq ~ 6.35m„d ~ 0.25tos. Since the volume 
and the quark mass are constant one can perform 
a continuum extrapolation from these points. In 
addition pion induced finite volume effects are un- 
der control as well, as they can be estimated by 
exp(— Lmps) ~ 0.03. 
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Table 1 

Estimates of parameters for the runs 

Results for our reference scale, the Sommer pa- 
rameter To, are presented in figure ^. 

For tuning the parameters to the values in ta- 
ble 1^ we further have to look at the pseudoscalar 
mass mps and the quark mass. This latter can 
be characterized by the dimensionlcss quantities 
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2 Gps 
and 

Mr = {r^mpsf. (6) 

The two definitions become proportional to each 
other for small masses. First results are shown 
in figures ^ and ^. In figure ^ there are two hor- 
izontal lines indicating TOs and the mass region 
we want to reach, ^to^- Depending on the size of 
the mass 8^ x 16, 12^ x 24 and 16^ x 32 lattices 
have been used. 
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Figure 2. Sommer scale tq for (3 = 5.1 



4. CONCLUSIONS AND OUTLOOK 

We presented preliminary results from simula- 
tions with three dynamical quarks. The results 
show that it should be possible to reach small 
quark masses within the proposed settings. For 
the moment we are exploring the region below the 
strange quark mass. 
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Figure 4. Two definitions of the quark mass 
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